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ABSTRACT 

We establish a controlled comparison between the properties of galactic stellar halos 
obtained with hydrodynamical simulations and with particle tagging. Tagging is a fast way 
to obtain stellar population dynamics: instead of tracking gas and star formation, it ‘paints’ 
stars directly onto a suitably dehned subset of dark matter particles in a cohisionless, dark- 
matter-only simulation.Our study shows that there are conditions under which particle tagging 
generates good hts to the hydrodynamical stellar density prohles of a central Milky-Way-like 
galaxy and its most prominent substructure. Phase-space diffusion processes are crucial to 
reshaping the distribution of stars in infahing spheroidal systems and hence the hnal stellar 
halo. We conclude that the success of any particular tagging scheme hinges on this diffusion 
being taken into account, at a minimum by making use of ‘live’ tagging schemes, in which 
particles are regularly tagged throughout the evolution of a galaxy. 

Key words: Galactic halo - simulations: SPH. 


1 INTRODUCTION 

Observations of the stellar halo around local group galaxies provide 
strong constraints on different models of galaxy formation (for his¬ 
torical context see, for instance, |Eggen et al.|1962| |Searle & Zinn| 
|1978| ). In today’s standard ACDM cosmology, large galaxies such 
as the Milky Way are predicted to form in part through the progres¬ 
sive mergers of smaller progenitor galaxies ([White & Rees|1978| 
IFrenk et al.|i^ [White & Frenk|1991||Kauffmann et al.|1993t , 

and the observed stellar halo is formed, in one way or another, from 
the remnants of this hierarchical assembly (|Helmi & White] 1999 1 
IHelmi et al.|l^|Abadi et al.|2006| l. 

While the current consensus is that the stellar halo forms pri¬ 
marily via accretion of stars from tidally stripped satellites (e.g. 
[Bell et al.|2008| ) there remain a number of open questions: for in¬ 
stance, does the halo form exclusively from accretion, or are some 
halo stars formed within the main galaxy (referred to as in situ stars) 
during dissipative collapse? Are further stars kicked up from the 
disc to form an important fraction of the halo ( [Zolotov et al.|200^ 
[2010[[Font et al.|20TT|[Tissera et al.|2013[ l? Do the main progeni¬ 
tors of the stellar halo survive as satellite galaxies to the present day 


([Cooper et al.|201()| ) or are they mostly fully disrupted ( [Bullock &| 
[Johnston[2005[ )? 


In principle, cosmological simulations of galaxy formation 
which self-consistently follow the evolution of a baryonic compo¬ 
nent can be used to predict the properties of galactic haloes and 
provide answers to these questions. In practice, however, this is 
challenging because large particle numbers are required to resolve 
any detail within the faint, diffuse stellar halo, which contains only 
a few percent of all the stars in the galaxy. Only recently has it 
become possible to start studying some halo properties in this self- 
consistent way (e.g. [Zolotov et al.[2009[ l. 


As an alternative, some authors have proposed making use 
of dark matter (DM)-only N-body simulations - which are signifi¬ 
cantly less computationally expensive than hydrodynamical simu¬ 
lations - where stellar populations are ‘painted’ onto DM particles 
to reproduce the collisionless assembly of the stellar halo ( Bullock[ 


[& Johnston|2005[ [Cooper et al.|2010[[Rashkov et al.|2012| |. Such 

methods have been used to make quantitative predictions of halo 
substructure and dy namics ([Cooper et al.|201l| [Helmi et al.|2011[ 

[Cooper et al.|2013| [Gomez et af|2013[ ), but the validity of the as¬ 

sumptions underpinning these DM-only models remains controver¬ 
sial (e.g. IBailin et al.|20T4t . 


So at present there are two techniques for investigating stel- 
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lar halo structure - hydrodynamic simulations and particle tagging 
- that have their own strengths and weaknesses. Recent work has 
demonstrated that significant discrepancies between the two ap¬ 
proaches can arise ( [Bailin et al.||2014| t. In this paper, we will in¬ 
vestigate a specific tagging scheme (that of [Cooper et al.| 2010) and 
by comparing its predictions to those of hydrodynamical simula¬ 
tions, aim to understand better how and why differences between 
predictions arise. While this does not immediately resolve the ques¬ 
tion of how to produce “correct” predictions for the stellar halo, it 
does give a physical basis for understanding discrepancies and so 
highlights some essential prerequisites for realistic modelling. 

This paper is structured as follows. In Section 2, we provide an 
overview of the tagging technique, stating its main assumptions; in 
Section 3, we describe the simulations used in our study, and how 
they are used to establish a controlled comparison between tagging 
and SPH; in Section 4, we present the outcome of the simulations, 
and compare the structure of the different stellar haloes obtained; 
in Section 5, we discuss the role of diffusion processes in shaping 
realistic tagged haloes, and give a new interpretation for the model- 
to-model discrepancies found in the literature. Finally, in Section 6, 
we explain the overall physical picture that emerges and highlight 
how this should affect the future direction of tagging and broader 
investigations of the stellar halo. 


ging cannot be accurate, especially considering also the limitation 
(ii). Secondly, baryonic discs can cause real haloes to be less pro¬ 
late than the tagged DM haloes ( [Bailin et~^|2014| ). Finally, the 
presence of a disc or of gas expulsion by SN feedback could also 
affect the dynamics of the DM component of a hydrodynamic sim¬ 


ulation - for instance creating cores in satellites ( [Pontzen & Gov-[ 
[ernato[20T4] ). Dynamical core creation affects the orbit of stars just 


as it does dark matter ([Navarro et al. 1996 

[Read & Gilmore|2005| 

[Pontzen & Governato|2012|[Maxwell et al. 

2012|) and, furthermore. 


the cored satellites are then more easily disrupted ( [Penarrubia et al.[ 
[2010||PS-ry et al. [20 1 2[ [Zolotov et al.|2012[ ). Conversely, baryonic 


dissipation could give rise to a contraction of the host subhalo mak¬ 
ing it more resilient to tidal disruption. 

The specific tagging scheme we investigate in this paper is de¬ 
scribed by [Cooper et al.[ ( [20T0] ), hereafter CIO, where a fixed most- 
bound fraction (/mb, usually chosen to lie between 1 and 10%) of 
DM particles in each halo are tagged at every simulation output 
time, being assigned ‘stellar masses’ (as well as ages and metal- 
licities) according to the prescriptions of a semi-analytic model of 
galaxy formation, GALFORM (see [Cole et al.|200()| for details and 
[Baugh|20'06| for an overview of hierarchical galaxy formation with 
semi-analytic models). 


2 PARTICLE TAGGING 


3 SIMULATIONS 


Various methods have been proposed for associating a stellar com¬ 
ponent with particles from D M-only simulations ([Bullock & John" 


ston[2005|[Cooper et al.[2010[[Libeskind et al.[201 1[ Rashkov et al. 

2012[ [Bailin et al.[[2014j ). DM particles that are tightly bound to 

their haloes are assigned ‘stellar masses’ according to an assumed 
star formation rate, and the evolution of these painted particles is 
then traced up to z = 0, where their final distribution is taken to 
represent the distribution of stars in the real stellar halo. 

These methods rely on three assumptions: 


(i) Stars form tightly bound to their parent halos (e.g. with bind¬ 
ing energies typically larger than 90 - 99% of the DM particles), 
with energy distributions similar to those of the ‘most-bound’ frac¬ 
tion of DM particles. 

(ii) Recently formed star particles and their tagged DM ana¬ 
logues subsequently follow similar phase space trajectories. In 
other words, selecting the correct initial binding energy for the 
tagged DM particles is assumed to be sufficient to ensure the cor¬ 
rect subsequent kinematics. This assumption could fail given that, 
for instance, DM particles and stars with similar binding energy 
will not necessarily have similar angular momentum. 

(iii) Baryonic effects are not important in shaping the stellar 
halo. For instance, the tidal disruption of subhalos is assumed to 
be unaffected by any interaction between baryons and dark matter 
and the presence of a baryonic disc is typically assumed not to af¬ 
fect significantly the distribution of stars in the halo (although some 
tagging schemes do include effects from an analytic disc potential, 
e.g. [Bullock & Johnston|2005] l. 

The first assumption is not likely to be problematic, but the 
second will fail if the stars are supposed to lie in a rotationally- 
supported disc; consequently only the accreted stellar halo is usu¬ 
ally considered in tagging schemes. The last assumption could af¬ 
fect the validity of tagging in several ways. Firstly, if in reality stars 
kicked up from the disc form a significant part of the stellar halo 
(see for instance [Zolotov et al.[[2009[ [Font et ar][201I[ |, then tag- 


We have made use of two sets of simulations, which we call the 
‘Durham’ and ‘Seattle’ simulations. Both sets contain smooth- 
particle-hydrodynamics (SPH) and DM-only (DM0) simulations 
of a Milky-Way like galaxy and its immediate environment. They 
use the volume renormalization or ‘zoom-in’ technique ( [Katz &[ 
[White[[T9^ to increase computational efficiency, with the mass 
and spatial resolution of the simulation decreasing with radius from 
the central galaxy. Each adopts a ACDM cosmology; the specific 
parameters adopted differ slightly with the Seattle (Durham) val¬ 
ues being = 0.24 (0.25), Da = 0.76 (0.75), Db = 0.042 
(0.045), as — 0.77 (0.9), ris — 0.96 (1.0) and in both cases 
Hq — 100/ikms“^ Mpc“^ = 73kms“^ Mpc“^. The setup of 
the two simulations is broadly comparable but crucial differences 
appear in the handling of star formation feedback as we will discuss 
below. 

The ‘Seattle’ galaxy, sometimes known as h258, was first run 
by [Zolotov et al.[ ( [2009] ) and then by [Zolotov et al.[ pOI2[ ) at the 
higher resolution employed in our current work. Its properties have 
been extensively discussed in the con text of the | Boylan-Kolchin[ 


[et al.[ ( [20TT] ) ‘too-big-to-faiT problem ( [Zolotov et al.|2012 


it has 

also been used as a case study of how baryonic discs can regrow 
following a major merger at z = 1 ( [Govemato et al.[[2009[ ). At 
the end of the simulation it has built a relatively low mass Milky 
Way (M 200 = 9 X 10^^M©, where M 200 is the mass enclosed in a 
sphere with 200 times the mean cosmic density) with a resolution 
of 1.3 X lO^M© (DM) and 2.7 x lO^M© (gas) and force softening 
of 170 physical parsecs. 

This simulation was run using GASOLINE ( [Wadsley et al.[ 
[2004[ ), a Tree-SPH code. When hydrodynamics is switched on, it 
uses a pressure-entropy interpolated kernel to eliminate artificial 
surface tension. Cooling, molecular hydrogen and star formation 
physics are described in [Christensen et aL]p012j ); feedback is im¬ 
plemented according to the [Stinson et al.[ ( |2006| ) blastwave model. 
Self-bound substructures are identified using Amiga’s Halo Finder 
( [Knollmann & Knebe|2009| l, hereafter AHF, although we processed 
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Figure 1. Projected stellar mass surface density at 2 ; = 0 for SPH stars and the ‘particle tagging’ outputs for the Seattle simulations. In these panels, tagging 
is always applied to the top five percent most bound particles (/mb= 5 %)- Note that a one-to-one correspondence between the DMO tagged and SPH tagged is 
not expected, as small differences between the SPH and DMO runs, combined with the stochastic nature of halo formation, will make them differ qualitatively 
(see Section!^. 


the final snapshot with SUBFIND (see below) to ensure the same 
structures are identified by both algorithms. 

The Durham simulation is based on an early version of the 
PM-Tree-SPH code Gadget-3, and takes the same inital condi¬ 
tions as halo Aq-C of the Aquarius suite ( |Springel et al.|2008] ). Its 
highest particle mass resolution (in a ~ 5/i“^Mpc region around 
the target halo) is similar to that of the ‘level 4’ simulation set in 
Aquarius. The final mass is M200 = 1.8 x with a reso¬ 

lution of 2.6 X 10^ M© (DM) and 5.8 x lO^M© (gas). The force 
softening is 260 physical parsecs. 

Durham baryonic processes are modelled as described in 
[Okamoto et al.| ( |2010| ), with a number of modifications designed 
to improve the treatment of supernovae-driven winds, which are 
explained in detail in |Parry et ^ { |2012| l. Additional details of 
this simulation are also presented in [Okamoto] ( |2013| ). The cen¬ 
tral galaxy contains a massive centrifugally-supported disc as well 
as a dispersion-supported spheroid. Self-bound substructures were 
identified with the SUBFIND algorithm ( [Springel et al.||2001| | as 
modified by |DoIag et"^ ( |2009| ) to take into account the internal 
energy of gas particles when computing particle binding energies. 

The central galaxy in the Durham simulation (of a ~ 
5/i“^Mpc region around the target halo) contains a massive 
centrifugally-supported disc as well as a dispersion-supported 
spheroid. The final mass of its halo is M200 = 1.8 x 10^^ M©. 
This simulation was used by [Parry et aL| ( |201^ to study the satel¬ 
lite system of a Milky Way-like system. It provided a good match 
to the average of the satellite luminosity functions of the MW and 
M31 but formed its brightest satellites in excessively massive ha¬ 
los, anticipating what became known as the too-big-to-fail problem 


[Boylan-Kolchin et ^ { [201 Ij t. One of the satellites in this simulation 
generated a core in its halo as a result of a large infiow of gas into 
the centre triggered by a triple merger and its subsequent violent 
expulsion in a starburst. 

As described above there are a number of numerical and phys¬ 
ical differences between the two sets of simulations. We will see 
below that the most important distinction is the type of stellar feed¬ 
back implemented. In the Durham case, feedback should be seen as 
‘passive’ in the sense of jPontzen & Governato| ( [20T4) ; it produces 
little or no coupling between the baryonic and DM components 
of the simulation. Conversely the Seattle simulations have ‘active’ 
feedback which demonstrably couples baryons and DM with en¬ 
ergy being passed from the former to the latter ( jPontzen & Cover- j 
|nato|2012] l. This means that, for instance, some of the main satel¬ 
lites in the Seattle simulations develop central density cores ( |ZoIo-| 
[toy et aL|2012j l, which is not generally the case in the Durham sim¬ 
ulation, with the exception of the single satellite mentioned above 
which, in any case, was disrupted shortly after being accreted into 
the main halo. Thus, when comparing haloes and tagged DM pop¬ 
ulations in these two sets, we can appreciate the effect of the very 
different feedback recipes implemented. 

3.1 Post-processing 

In this paper, we will establish a controlled comparison between the 
self-consistently formed stars in the simulations described above 
and “tagged” stars in dark-matter-only (DMO) runs of the same 
regions. For full insight we will also need to study stellar haloes 
obtained through algorithms intermediate between these two cases. 
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Figure 2. Spherically averaged density profile of the main halo, excluding self-bound substructures, for the SPH stars and tagged DM, for most bound fractions 
(/mfos) of one and five per cent, in the Seattle (left panel) and Durham (right panel) simulations. The five per cent cases produce output densities within a 
factor of a few of the SPH results over seven orders of magnitude in stellar density in all cases. 


Our series of approaches are similar to those used by [Bailin et^ 
( |2014| ), but differ fundamentally in that they all employ ‘live’ star 
formation - associating stellar mass with dark matter particles in a 
time-dependent way. They are labelled as follows: 

• SPH STARS - the stellar component of the SPH simulations. 
These represent an obvious basis from which to assess the success 
of a tagging scheme, although one should bear in mind that there 
are a number of differences between SPH and real galaxies; see for 
instance [Stinson et al.| ( |2010| ) ; [Creasey et al.| ( |2011| ) . 

• SPH TAGGED SF-REAL - DM particles are tagged within 
the SPH simulation, and the stellar masses assigned to the parti¬ 
cles are calculated using the actual star formation rates (SFR) of 
the SPH galaxies. Thus the star particles are not used explicitly but 
SFRs are guaranteed to match and baryonic effects such as feed¬ 
back and the presence of a disc can still affect the tagged DM’s 
dynamics. 

• SPH TAGGED SF-MODEL - as above, DM particles are 
tagged in the SPH simulation; but SFRs are calculated using an 
analytic prescription. Specifically, we assume a power law relation 
between halo DM mass and SFR. This relation is allowed to vary 
with redshift, and the power law indices are obtained by fitting to 
the SPH SFR as a function of halo mass. Note that CIO’s tagging 
scheme actually uses the semi-analytic model GALFORM to ob¬ 
tain SFRs, but in this paper we adopt a simpler prescription to ob¬ 
tain SFRs consistent with, but no longer identical to, those of the 
SPH comparison simulations. 

• DM0 TAGGED SF-MODEL - particles are tagged in a DM- 
only (DM0) run of the galaxy. Ultimately, it is the validity of 
this approach that we wish to investigate, since it resembles most 
closely the scheme used by CIO to produce observable predictions 


of halo properties. The SFRs are calculated in the same way as the 
SPH TAGGED SF-MODEL run. Note that again we do not use the 
CIO model to obtain the SFRs since we wish to test the physics of 
tagging, not of semi-analytic galaxy formation models. 

If tagging fails due to an intrinsic difference in the kinemat¬ 
ics of DM and star particles, then this should be clear when com¬ 
paring the haloes obtained in the SPH STARS and SPH TAGGED 
SF-REAL. If, on the other hand, it fails because baryonic effects 
do affect significantly the assembly of the stellar halo, then the 
disagreements should arise when comparing the DM0 TAGGED 
SF-MODEL to the SPH STARS and SPH TAGGED SF-REAL 
runs. Finally, the SPH TAGGED SF-MODEL is a control run to 
check that any disagreements between the SPH and DM-only tag¬ 
ging do not in fact occur because the scaling relation used to pre¬ 
scribe the star formation rates is too naive. 


4 HALO STRUCTURE 

The principal aim of particle tagging is to produce tagged stel¬ 
lar haloes with realistic substructures at z = 0 (such as stellar 
streams) and having densities and dynamics comparable to those 
of real haloes. The purpose of the present paper is to investigate 
some of the physics relevant to stellar halo formation rather than 
to tackle observational matters - this will be left to a companion 
paper. Cooper et al (in prep) - but we nonetheless need to start with 
a brief comparison of the observable outputs. 

To obtain a qualitative view of how well DM tagging repro¬ 
duces the features of the SPH stellar halo, we obtained 2-D pro¬ 
jected density maps of the SPH stars and their tagged DM ana¬ 
logues at z = 0, tagging first the DM in the SPH, and then in the 
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Figure 3. Spherically averaged density profiles of the four most massive satellites for the SPH stars and the different tagging realizations. Densities are 
normalized so that the total tagged mass equals the SPH stellar mass. Radii are normalized to rmax, the radius at which the circular velocity achieves 
its maximum value. The bottom right panel shows the profile of the largest satellite in the Durham simulations. All other satellites come from the Seattle 
simulations. For Seattle Satellite 1, the scale lengths are significantly underestimated in all DMO cases because the feedback coupling - that causes collisionless 
particles including stars to migrate outwards - is missing. In the other cases reasonable agreement can be achieved using 5% tagging. 


DM-only simulation (Figure 0. We find that the SPH TAGGED 
outputs produce a stellar halo containing much the same substruc¬ 
ture as the SPH stars: in both cases, stellar streams are clearly 
discernible, associated with tidally stripped satellites of the main 
galaxy, and seem to have comparable densities. This suggests that 
the SPH tagged DM and SPH stars have similar density distribu¬ 
tions by the time the simulations reach z = 0 - though this does 
not necessarily imply that the SPH stars and tagged DM were in 
agreement at the time of tagging. This is an important distinction 
which will be discussed further in the next section. 

The resolution in the tagged DM tends to be lower than in the 
SPH stars. This is not an intrinsic limitation of particle tagging, but 
an artefact of the way we carry out the comparison: here we tag 


DM in an SPH simulation (and a corresponding DM-only simula¬ 
tion), where the mass and spatial resolution in the DM component 
is lower by a factor of about 10 than those achieved by recent pure 
A^-body codes (such as the Aquarius haloes, on which CIO applied 
their tagging scheme). In order to compare the tagged and SPH star 
images more fairly, one can thus downsample the latter, taking only 
every fifth particle (top right panel of Figure [^. The effect of this 
is to decrease the resolution of the streams, which is also observed 
in the SPH TAGGED. In other words, the fact that substructure is 
somewhat ‘blurred’ in the tagging images is related to the artifi¬ 
cially low resolution. 

In fact the resolution issues become even more subtle and can 
be seen to be responsible for the visible tendency of tagging to pro- 
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duce larger abundances of subhaloes, even in the SPH TAGGED 
SF-REAL case where the effects of baryonic feedback on SFRs 
are directly taken into account, and where tidal interactions with 
the galactic disc still have the ability to disrupt satellite haloes. The 
number of dark matter particles tagged is independent of the star 
formation rate; therefore, in a halo with low SFR, a large number 
of particles may be tagged with very small stellar masses. In this sit¬ 
uation the mass resolution locally becomes much higher than that 
of the SPH stars, which have fixed mass per particle at birth irre¬ 
spective of SFR. So, although overall the resolution is lower in the 
tagged halo, we verified that for small subhaloes the resolution is 
in fact larger in the tagging realizations than in the SPH stars. 

There seems to be fairly little qualitative difference between 
the SPH TAGGED SF-REAL and SPH TAGGED SF-MODEL 


runs, which suggests that the simple prescription used for assigning 
SFRs to DM haloes is adequate for our purposes. In the case of the 
DM0 tagged simulation a qualitative comparison is less straight¬ 
forward, as small differences between the SPH and the DM0 sim¬ 
ulation runs, combined with the chaotic nature of halo structure for¬ 
mation, means that even the locations of DM substructures are quite 
different. Although streams and substructures are present, they do 
not correspond directly to the SPH streams, and there seem to be a 
significantly larger number of subhaloes than in the SPH case. We 
will argue below that this last difference may, at least in some cases, 
be an effect of the ‘bursty’ feedback implemented in the Seattle 
simulations, where some satellites develop DM cores in their in¬ 
ner regions, and thus become more susceptible to tidal disruption 
when passing through the galactic disc - see [Zolotov et al.| ( [2012| ); 
[Brooks & Zolotov| ( |2014| ) and also [Penarrubia et al.[ ( [2010[ ). 

For a more quantitative comparison between the SPH stars and 
tagged haloes, we turn to spherically-averaged density profiles for 
the main halo and a few satellite galaxies. Figure]^ shows the main 
stellar halo density for various tagging realizations for the SPH and 
DM0 cases (top and bottom panels) in the Seattle and Durham sim¬ 
ulations (left and right). In each panel, the thick black dashed line 
shows the SPH STARS. For all three tagging realizations, taking 
fmb = 5% produces output densities within a factor of a few of 
the SPH results over seven orders of magnitude in stellar density. 
Notably, the Seattle SPH TAGGED cases are in particularly good 
agreement - this result can be tied to strong active star formation 
feedback and will be discussed further below and in SectionjS] 

Conversely the worst agreement arises with DM0 tagging 
within the inner few kpc of the main galaxy. Tagging will be a 
poor description of the star distribution in regions where baryons 
dominate, i.e. near the disc and bulge of the main galaxy. The in¬ 
accuracies can then propagate outwards if a significant fraction of 
the halo stars were originally disc stars that were scattered out to 
the stellar halo. A more detailed discussion of this process and its 
importance in Seattle simulations can be found in [Zolotov et al.[ 
( [ 250 ^ . 


We also obtained density profiles for the main satellites (Fig¬ 
ure 1^. In reading order, the figure shows the first to third most 
massive satellite in the Seattle simulations, then the most massive 
of the Durham satellites. Once again the SPH stars are shown by 
the thick dashed line - one can read off that the agreement between 
the SPH stars and tagging fmb = 5% in both SPH and DM0 is 
broadly quite good. However the most massive Seattle satellite (top 
left panel) shows a striking discrepancy between the DM0 and SPH 
cases. It has previously been shown for the Seattle simulations that 
the SN feedback causes the dynamics of the DM and baryonic com¬ 
ponent to couple strongly, creating DM density cores ( [Pontzen ^ 
[Governato|2012[ | in the central regions. Tagged and ‘real’ star par- 
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Figure 4. Evolution of the half-mass radius of different tagged populations 
and the stars they are meant to represent in the Seattle main galaxy (top 
panel) and its three largest satellites (next 3 panels, in order of satellite size). 
The lowest panel shows Durham satellite 1. The colour scheme is as for the 
density profiles: black for the SPH stars, blue and red for the SPH tagged 1 
and tagged 5 respectively, and cyan and magenta for the DM-only tagged 1 
and tagged 5 respectively. In all cases, the half-mass radius at early times is 
sensitive to the tagging mode, but at late times this dependence is substan¬ 
tially reduced. This is tied to diffusion of particle energies as described in 
the text, and is a particularly strong effect when ‘active’ feedback is present 
as in the Seattle SPH simulations. 

tides alike will be thrown to larger radii by this process. But in 
the DM0 simulation, the process is of course missing - so the DM0 
TAGGED cases form stellar components with significantly too short 
a scale length. This begins to show that the agreement between tag¬ 
ging and SPH will depend on what kind of feedback recipe is used 
in the SPH (i.e. ‘active’ vs. ‘passive’). 

In fact, this deep connection between the mode of stellar feed¬ 
back and the behaviour of tagging algorithms can be further re¬ 
inforced by studying the difference between the different tagging 
fractions (fmb =1% or 5%, hereafter “tag 1” and “tag 5”). As 
we commented briefly above, in the Seattle SPH TAGGED out¬ 
puts setting 1% or 5% seems to make surprisingly little difference 
to the final density profiles of the halo, especially in the satellites, 
or in the main galaxy in the outer regions. Increasing the number 
of tagged particles by a factor of 5 does not for instance change 
the scale length of the exponential density profiles, or their general 
shape. This suggests that the overall density profiles are primarily 
being set by dynamical processes after tagging, not by the details 
of the tagging itself. The expected increased prominence of these 
processes when feedback is active is reflected in near-complete in¬ 
sensitivity to fmb in the Seattle SPH TAGGED cases. 

The Durham simulations have a more passive type of feed¬ 
back and are therefore missing the dynamic redistribution of stars 
that we are claiming. The bottom right panel of Figurej^shows the 
profile of the largest satellite in the Durham simulations: agreement 
between the tagging realizations and the SPH stars is not as good 
as for the Seattle satellites. Moreover, the SPH tag 1 and tag 5 runs 
do not resemble one another as closely as they do in the Seattle 
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Figure 5. Diffusion in the rank binding energy distribution of recently formed stars and their analogue tagged DM in the Seattle main halo (left panel) and 
largest satellite (middle panel), as well as the Durham largest satellite (right panel). From top to bottom, the panels represent the simulations at simulation 
output times: t = 800 Myr, t = 3200, t = 4500, where the stars are those formed at t = 800 Myr, and the DM has been assigned a stellar mass at the same time. 
The evolution over time is characterised by a drift to larger scale radii, and (in most cases) towards closer agreement between initially distinct populations. As 
also seen in Figure[^ the effect is strongest in the Seattle SPH cases where active feedback generates strong potential fluctuations. 


satellites. This is suggestive, but we now turn to a more explicit 
way to verify that diffusion of stars through the energy space is a 
significant effect. 


5 DIFFUSION 

In the previous section we compared the stellar haloes from our 
various tagging realizations and found, from a number of perspec¬ 
tives, agreement between apparently very different approaches. We 
have suggested that these can be tied to dynamical redistribution of 
the stars after they have been formed; we will now provide a direct 
analysis to support that claim. 

We first trace the evolution of the radial scale length of a pop¬ 
ulation of tagged particles over a period of a few Gyr after its initial 
assignment, as follows. In an early simulation snapshot {z = 8), we 
select the most-bound 1 and 5 percent of DM particles, as well as 
the recently formed stars. We then obtain the radius enclosing half 
of the total mass of these particles. We track these same particles 
over time, calculating their half-mass radius at each snapshot. The 
results are shown in Figure for (top to bottom panels) the ma¬ 
jor progenitor and the three most massive satellites in the Seattle 
simulations, followed by the most massive satellite in the Durham 
simulations. 

At the tagging snapshot, z = 8 (far left of each panel), the 
half-mass radii do not typically match those of the stars, with initial 
discrepancies between SPH TAGGED SF-MODEL (dotted lines), 
DM TAGGED SF-MODEL (solid lines) and SPH STARS (dashed 
line) reaching an order of magnitude. The difference between tag¬ 
ging the 1 or 5 percent most-bound DM particles in the halos is 
also clearly discernable: as naively expected, when a larger number 
(by a factor of 5) of particles are tagged, their distribution is more 
diffuse, and the half-mass radius is larger. 

However, as we follow the evolution of this distribution, all 
populations “diffuse” outwards to larger radii; initially compact 
populations tend to diffuse faster and so the initial differences in 


scale radii are significantly eroded. Particularly in the Seattle cases, 
the tagged particles of the SPH runs (dotted lines) tend to converge 
to mimic the stars (dashed line) - irrespective of the binding ener¬ 
gies at which they were originally selected. In these cases the con¬ 
vergence between the half-mass radii of the tag 1, tag 5 and stars 
occurs on a timescale of around 3 Gyr, and these three initially dis¬ 
tinct populations become virtually indistinguishable. 

In the DM-only tagged assignments (solid lines) convergence 
is slower, but still significant. In Seattle satellites 2 and 3, for ex¬ 
ample, the half-mass radius discrepancy narrows from an initial dis¬ 
crepancy of a factor 2 at z = 8 to agreement by z = 0. However, 
unlike in the SPH case, it takes the entire Hubble time for this to oc¬ 
cur and in the case of satellite 2 it does not even run to completion. 
Furthermore in the central galaxy the diffusion process ceases after 
a major merger at a time of around 4 Gyr. This is consistent with 
a picture of diffusion being driven primarily by time-dependence 
in the potential of a spheroidal halo - after this time the galaxy is 
quiescent and, in the case of the SPH run, has developed a stable 
disc structure. 

For each satellite, the diffusion is faster and more sustained in 
the SPH case than in the DM0 case. This follows because baryonic 
feedback enhances the process by producing more stochastic time- 
dependence in the potential, particularly in the case of satellite 1 
(second panel from top in Figure]^ where the populations from the 
SPH simulation continue migrating outwards throughout cosmic 
time. As discussed in Sectionthis satellite produces a feedback- 
driven dark matter density core in the hydrodynamic cases. The 
outward migration is a symptom of this process and by definition 
is not present in the DM0 cases. 

With the above in mind, we can start to see the diffusion pro¬ 
cess as a double-edged sword. On the one hand it does tend to drive 
tagging realizations into better agreement, because it reduces sen¬ 
sitivity to the initial choice of which particles to tag. On the other, 
the strength of the diffusion can be enhanced by baryonic processes 
leading to a new, fundamental source of discrepancy with DM0 
simulations. In the Durham cases (illustrated here by Satellite 1, 
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the lowest panel of Figure]^ the DM0 and SPH cases are in much 
better agreement, but the convergence is also much weaker. Both 
these aspects follow from the relatively passive approach to feed¬ 
back taken. 

In practice, this means that - if one believes in the kind of 
bursty star formation histories that lead to the ‘active’ feedback 
coupling - DM0 tagging schemes will need to tag particles that are 
less bound than the regions in which star formation actually takes 
place, because they miss processes that kick stars up to higher en¬ 
ergies. Concretely, the lower left panel of Figurej^ shows 5 percent 
tagging in the Seattle DM0 simulations produce better fits to the 
SPH star density profiles than 1 percent tagging - and yet Figure 
1^ shows directly that at the time of formation, the SPH STARS are 
more tightly bound even than the 1 percent case. 

5.1 Energy and angular momentum 

We can obtain another view of the importance of diffusion by 
studying the process in energy space. First, we obtain the energy 
distributions of recently formed stars and their tagged DM ana¬ 
logues (using different values of /mb) at the time the DM is as¬ 
signed a stellar mass, and track the evolution of these distributions 
(for the same set of particles) over a few gigayears. We confirm 
that, at the time of tagging, the recent stars and most-bound DM 
particles do have very different energies (top panels. Figure the 
three panels from left to right show Seattle MW, Satellite 1 and 
Durham Satellite 1). 

One can now see explicitly how in the Seattle SPH cases, 
the baryonic feedback drives the full SPH STARS and SPH TAG 
MODEL-SF energy distributions into agreement, irrespective of the 
starting /mb- Once again the DM0 cases do evolve to reduce ini¬ 
tially stark differences in the energy distributions, but this diffusion 
does not complete and the energy distributions of this population 
of stars are significantly different from those of the SPH case. 

Comparing satellite 1 from the Seattle and Durham SPH simu¬ 
lations (middle and right panel respectively), one can see explicitly 
how the difference between active and passive feedback is crucial. 
In the Durham case, the differences in the diffusion between SPH 
and DM0 are much less marked because the gas evolution has rel¬ 
atively little impact on the collisionless dynamics. Once again we 
see the dual impact of diffusion - it substantially reduces sensitivity 
to the initial distribution of star formation events, but the extent of 
diffusion depends on the nature of feedback. 

Finally we verified that a very similar view can be obtained 
in angular momentum space (Figure]^ - the chaos associated with 
potential fluctuations from mergers or stellar feedback causes the 
distributions to drift away from their initial arbitrary shapes. Only 
when stars form in a stable disc is the diffusion in both angular 
momentum and energy stopped. 

In conclusion, there is very strong evidence for the diffusion 
of orbits through phase space being a crucial ingredient in deciding 
the final distribution of stars in spheroidal systems. This in turn has 
significant implications for interpreting the existing literature and 
for future schemes attempting particle tagging, as we now discuss. 


6 CONCLUSIONS AND DISCUSSION 

We have re-implemented and investigated a particle tagging 
scheme for simulating stellar haloes in DM-only simulations, first 
used and described in [Cooper et al.| ( [20T0] t. The scheme assigns a 
‘stellar mass’ to the most-bound DM particles in each halo at every 
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Figure 6. The angular momentum distribution of a population of stars and 
their analogue tagged DM in the main halo, selected at high redshift (z = 8) 
in the Seattle simulations. Stars like these that end up in the halo, rather than 
the disk, have distributions of angular momentum evolves rapidly - similar 
to how their energy distribution changes (Figure]^. Stars forming later and 
contributing to the disk in the SPH case retain their angular momentum. 

simulation output time. In our case, we use this method to establish 
a comparison similar in spirit to that of jBailin et al.| ( [20T4] ), con¬ 
structing a series of different stellar halos which increasingly differ 
from the SPH stars. Our focus is on understanding the physics that 
sets the distribution of stars within each subhalo, and hence in the 
stellar halo of a Milky Way-like galaxy. 

Our main conclusions are: 

• Tagging a fixed fraction of dark matter within SPH simula¬ 
tions reproduces qualitative features of the SPH stellar halo struc¬ 
ture and sub-structure (Figure [^; tagged stellar haloes have den¬ 
sity profiles very similar to the SPH haloes (compare red and blue 
’tagged’ lines to the black dashed line in Figure]^. 

• Performing the same experiment in the DM0 run produces 
more pronounced differences between the tagged halo and refer¬ 
ence SPH stars halo, although for 5% tagging one obtains similar 
scale-lengths and density profiles that deviate from the SPH stars 
by less than one dex (see Figurej^. The exception to this pattern is 
Seattle’s satellite 1, which has a large dark matter core in the SPH 
run; therefore, by definition, collisionless particles in the DM0 run 
end up on the ‘wrong’ orbits relative to the SPH case. 

• All the above results can be tied to diffusion of the tagged DM 
particles through energy space. Populations which are tagged in 
differing regions nonetheless end up with similar phase space dis¬ 
tributions if sufficient time is given between the initial assignment 
and the point at which the stellar halo properties (densities, satellite 
profiles and so on) are measured (e.g. Figures]^ and [^. 

• The agreement between tagging and SPH is affected by the 
type of SN feedback recipes implemented in the SPH sub-grid 
physics. The Seattle simulations have a strong, ‘bursty’ feedback 
recipe ( [Stinson et al.||2006| ). We have already commented on the 
core in satellite 1, but even elsewhere in the Seattle simulations we 
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see considerably accelerated diffusion compared to the DM0 runs 
(e.g. Figure [^. Conversely the Durham simulations (which have 
‘passive’ feedback, rather than Seattle’s ‘active’ feedback; |Pontzen| 
|& Governato| 2014) drive the diffusion at rates much more similar 
to the DM0 case. 

Another way to phrase our conclusions is that (in a spheroidal, 
pres sure-supported structure) once a star has formed, the collision¬ 
less dynamics will make initially different phase space distributions 
of the tagged and star particles converge. Future work on stellar 
tagging needs to take this into account. Key recommendations are 
that (a) tagging must be done on-the-fly rather than at infall (see 
further discussion below); (b) one should consider tagging DM in 
regions considerably larger than the strict star formation region of 
a halo, to make up for missing energy kicks that might arise from 
astrophysics in the real universe. 

The recent work of |Bailin et al.| ( |201^ focused on comparing 
observable results between SPH stars and a somewhat different par¬ 
ticle tagging algorithm. The key difference between their work and 
ours is that their stars in a given satellite are tagged at the time of 
maximum mass (essentially just prior to infall), rather than dynam¬ 
ically throughout the simulation. The diffusion process therefore 
operates for a much shorter time in |Bailin et al.| ( |2014] ) and, accord¬ 
ingly, they see more stark differences between tagged and dynamic 
stars than we do. 

Diffusion itself may be driven in the DM-only simulation by 
mergers or tidal interactions ( [Stickley & Canalizo|2012[ |Kandrup| 
|et al.|200^ and accelerated in the SPH simulation by baryonic pro¬ 
cesses such as interaction with a disc or SN feedback ( [Valluri et al.| 
|2013| ). Further study of these processes will be necessary to un¬ 
derstand more fully the exact domain of applicability of particle 
tagging. 

Other proposals for improving and understanding tagging 
were given by [Bailin et"^ ( |2014| ) and focus on selecting a more 
appropriate region of the phase space, in particular by consider¬ 
ing angular momentum. Our work suggests that, in a spheroidal 
system, differences in initial angular momentum get smoothed out 
in much the same way as the differences in energy. Therefore se¬ 
lecting the ‘correct’ particles to tag is as much about studying dy¬ 
namical evolution as it is about characterising the orbits of recently 
formed stars. 
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